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ABSTRACT 

In this paper the transient thermal stress problem for a hollow 
elastic cylinder containing an internal circumferential edge crack is 
considered. It is assuned that the problem is axisymmetric with regard 
to the crack geometry and the loading, and that the inertia effects are 
negligible. The problem is solved for a cylinder which is suddenly cooled 
from inside. First the transient temperature and stress distributions 
in an uncracked cylinder are calculated. By using the equal and opposite 
of this thermal stress as the crack surface traction in the isothermal 
cylinder the crack problem is then solved and the stress intensity factor 
is calculated. The numerical results are obtained as a function of the 
Fourier nunber tD/b^ representing the time for various inner-to-outer 
radius ratios and relative crack depths, where D and b are respectively 
the coefficient of diffusivity and the outer radius of the cylinder. 


INTRODUCTION 

Cracking of brittle solids due to thermal stresses is a well-known 
phenomenon. In the absence of additional external loads, under thermal 
stresses, because of the self-equilibrating nature of the stress state, 
the cracking may not always lead to a through-thickness or catastrophic 
fracture. For example, in [1] it was shown that in a suddenly cooled 
hollow glass cylinder an axial initial flaw penetrated into the cylinder 
wall only partially but propagated axially the entire length of the 
cylinder as a part- through crack. Similarly, one would expect that if 
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the Initial flaw were to be circunferential a part-through crack could form 
and propagate along the entire circumference of the cylinder. In either 
case, the depth of the part-through crack may be estimated by comparing 
crack arrest fracture toughness of the material with the stress intensity 
factor. This, in turn, requires the calculation of the stress intensity 
factor for the canbination of given mechanical, thermal, and residual 
stresses as a function of the crack depth. If the consideration is 
restricted to linearly elastic materials, then for each loading condi- 
tion the crack problem may be treated separately. 

In this paper only the transient thermal stress problem is considered. 
The problem is that of a long hollow circular cylinder containing an 
internal axisymmetric circunferential edge crack which is suddenly cooled 
from inside (Fig. 1, c=a, d<b). It is assumed that the resulting tran- 
sient thermal stress problem is quasi-static; that is, the inertia 
effects are negligible. Previous studies on dynamic thermoelasticity 
seem to bear out the validity of this assumption which, of course, 
simplify the problem quite considerably (see, for example, [2] and [3]). 
Also, all thermoelastic coupling effects and the possible temperature 
dependence of the thermoelastic constants are neglected. By taking advan- 
tage of the linearity of the material the thermal stress problem in the 
cracked cylinder is considered in two parts. The first problem is the 
evaluation of transient thermal stresses in a hollow cylinder without the 
crack. The second is the isothermal perturbation problem for the cracked 
cylinder in which the crack surface tractions equal and opposite to the 
thermal stresses obtained from the first problem are the only external 
loads. The superposition of the two solutions gives results for the 
thermal stress problem for the cracked cylinder. Needless to say, the 
important information with regard to fracture initiation and propagation 
in the cylinder is contained in the second problsn. 


THE THERMAL STRESSES 

In the solution of the basic crack problem given in [4] it was 
assumed that the z=0 plane is a plane of symmetry with respect to the 
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external loads as well as the crack geometry. Consequently, the crack 
surface traction was the only external load and the crack problan 

was one of Mode I. On the other hand, if z=0 is not a plane of symmetry, 
then the shear tractions and on the crack surfaces would not be 
zero, the problem would become one of mixed mode, and extremely compli- 
cated. Even though any quasi-static transient thermal stress problem 
giving 

^^^(r, 9 , 0, t ) — f(r,9,t), ~ Tg 2 (r, 0 ,O,t) — 0 (1) 

for the uncracked cylinder can be solved by using the technique developed 
in [4] (with f as a known arbitrary function), in this paper, for sim- 
plicity, it is assumed that the thermal stresses are independent of 0 
and z. The thermal stress problem is, therefore, very simple (see, for 
example, [5], [6], [7]). Thus, if Tq is the initial temperature of the 
cylinder corresponding to zero stress state and T(r,t) is the temperature 
distribution at time t, then defining 

0(r,t) = T(r,t) - Tq , (2) 


the stress component of primary interest may be expressed as 

b 

^ 6(r,t)r dr - 0{r,t)] , 


(3) 


where E, v, and a are respectively the Young's modulus, the Poisson's 
ratio, and the coefficient of thermal expansion of the material. 

The temperature distribution is obtained by solving the diffusion 
equation 


V^0 


i li 

D 3t * 


(4) 


under the initial condition 9(r,0) = 0 and appropriate boundary condi- 
tions, where D is the coefficient of diffusion (i.e., D = k/pc; k, p, c 
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being the coefficient of heat conduction, the mass density and the 
specific heat). A simple set of boundary conditions which would lead to 
conservative results for many transient cooling problens may be stated as 

8(a,t) = (VTo)H(t) = eji(t) , (t>0) , (5) 


^e{b,t) = 0 , (t>0) 


( 6 ) 


where H(t) is the Heaviside function and for cooling problems e„ is a 
negative constant. 

Equation (4)may be solved bydefining 0 = e„ + u(r,t) and by using the 
standard technique of separation of variables. However, in this case, for 
small values of time (which is the period of main practical interest) 
the resulting series converges very slowly. An alternative technique 
much more suitable for our purpose is the use of Laplace transform. In 
this case, adequatley accurate representation of the solution for 
small times may be obtained by approximating the Laplace transform of 
e(r,t) by its asymptotic expansion for large values of the transform 
variable. The temperature distribution may thus be obtained as 


e« 


X ^ I*" 


-x^/4Dt 




-^;erTCL 


’2(Dt) 


where 


Li’s 
nVO^ 


-X„(dTT) 


-x_2/4Dt 
e " ] 


x-j = r-a, X 2 = 2b-r-a, X 2 = 2b+r-3a, x^ = 4b-4-3a, 


(7) 

( 8 ) 
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and the functions E^Cr), given by 

E-j “ Eg = -E 3 = “E^ = A^/B-j , 

F, = (A 2 B^-A^Bg)/B^^ , Fg = -(A^Bg+AgB^)/B ^2 ^ 

F 3 = (3A^Bg-AgB^)/B^ , F^ = (3A^Bg+AgB,)/B,^ , 

G] = (B^A3-A^B3-AgBg)/B^2, Gg = {B^A 3 -A^B 3 +AgBg)/B^^ . 

G 3 = (-B^A3+A^B3+3AgBg)/B^^ - 2 A^Bg^/B^^ , 

64 = (-B^A3+A^B3-3AgBg)/B^^ - 2 A^Bg^/B^^ , 

A^ = (br)"^ , Ag = -(3r+b)/[8(br)^'^2] , 

A 3 = 3(5rf3b)(b-r)/[128(br)^/2] , 

B^ = (ba)"^ , Bg = -(3a+b)/[8(ba)^2] , 

B 3 = 3(5a+3b)(b-a)/[128(ba)^2] . ( 9 ) 


Substituting from (7) into (3) the axial thermal stress may then 
be calculated. This step is carried out numerically as the evaluation of 
the related integrals in closed form is rather difficult. 


THE CRACK PROBLEM 

The crack problem may be solved by using equal and opposite of the 
axial stress obtained from (3) as the crack surface traction in the 
cylinder containing a circumferential crack and by treating the problem 
as being isothermal and quasi-static. The details of the formulation 
of the general nonaxi symmetric problen is given in [4]. In the solution 
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given in [4] the r, e, and z components of the displacement vector are 
expressed in terms of the sum of a harmonic potential associated with an 
infinite elastic space containing a plane of symnetry and a set of four 
harmonic functions which are equivalent to Papkovich-Neuber potentials 
in cylindrical coordinates [8]. Fourier and Hankel transforms are used 
to formulate the problem. The elasticity problem (which is considered 
only for 0£Z«») is subject to the following boundary conditions 


o^^(a,z,t) = 0, a^^(a,z,t) = 0, (0£Z<», 0£t«») , (10) 
a^j,(b,z,t) = 0, a^^(b,z,t) = 0, (0<Z<», 0<t<») , (11) 
arzCi'iOjt) = 0* (air_<b, 0<t«»), (12) 
cTzz(>^»0»t) = -ajz (r,t) , (a<r<d, 0<t<») , 

Uz(i*»0»t) = 0, (d^r^b, 0<t<~) . (13) 


The homogeneous conditions (10)-(12) are used to eliminate some of the 
unknown functions arising from the integration of the differential equa- 
tions and the mixed boundary conditions (13) gives an integral equation. 
By defining the derivative of the crack surface displacement as the 
unknown function, namely 

Uz(r,0‘‘‘,t) = (Ji(r,t) , (a<r<d, 0£t<») , (14) 

after a rather lengthy analysis the integral equation for ^ was found to 
be [4] 

a 

+ l2(s,r) - L2(a,r)}<}.(s,t)ds = - ^ <^zz(r,t) , 

(a<r<d, 0<t<«) , (15) 
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where p is the shear modulus and the kernels L-j and Lg are given in [4]. 
Note that in this problon time t enters into the analysis through 
only and (15) must be solved for each value of t separately. The integral 
equation (15) is singular and corresponds to an "edge crack" problem. Its 
nimerical solution may be obtained by using a Gaussian integration pro- 
cedure in a relatively straightforward manner^ ' [9]. The solution of 
(15) is of the following form [10]: 

4 »(r,t) = , (a<r<d) , (16) 

where f is a bounded function. 

The quantity of primary interest here is in the stress intensity 
factor k which is the fracture mechanics parameter and is defined by 

k = lim >^^(r-d) a„(r,0,t) . (17) 

r^ 

After solving the integral equation k may be obtained from 

k = lim v^Cd-r) (-j^) <j) (r,t) . ' (18) 

r-Kl ' 

RESULTS 

The numerical results giving the temperature e(r,t), the thermal 
stress a 22 (r,t), and the stress intensity factor k(d,t) have been obtained 
for the radius ratios a/b = 0.3, 0.5, 0.7, and 0.9 (Fig. 1), Time is 
represented through the dimensionless Fourier number which is defined by 

Fq = Dt/b2 (19) (*) 

(*)See [lO] for the treatment of an edge crack problem. 
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itfiere D, t, and b are the coefficient of diffusivity, time, and the 
outer radius of the cylinder, respectively. 

The calculated tanperature distribution and the corresponding axial 
stress in the cylinder are shown in Figures 2-9. Note that in the cool- 
ing problem the normalizing temperature 0 = T„-T and 0 = T-T. are both 
negative. As time goes to infinity the temperature ratio 0 / 0 ® approaches 
one and the stress approaches zero. The figures show that these limits 
are approached much faster for the larger values of a/b, that is, for 
relatively thinner cylinders. The peak value of the stress which occurs 
at t=0 has the expected value of Ea0»/(l-v). This is the thermal stress 
in a fully constrained plate undergoing a sudden uniform temperature 
change 0«. 

The stress intensity factor obtained from (15) and (18) and normalized 
as 

k (d) = — — . (20) 

‘ Ea0„/e/(l-v) 

is shown in Tables 1-5, where I = d-a is the crack depth (Fig, 1). Note 
that the thermal stress is statically self-equilibrating. Therefore, 
for £<h = b-a the resultant force on the crack surface is compressive 
and its magnitude decreases with increasing 1. Consequently, at a given 
time t (or for constant Fq) the stress intensity factor ratio kj will 
decrease as I increases. From the tables it may also be observed that 
for a given crack depth Z generally k^ first increases, goes through a 
maximum, and then decreases as F^ increases. This may be explained by 
the change in stress profile with increasing time. 
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Table T. Stress intensity factors for internal edge cracks 
subjected to transient thennal stresses Ckj(d) = 


/I 


l-v 

Ea0o 


= 0.1. Fo= 


Dt 




z 

h 

Ugijjjljll 


Fb = .001 

fb=.odr 


fb = .05 

kT(d) 

kj(d) 

kj(d) 

kjCd) 

1^ 

kjCd) 
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iSB 
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gm 
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BH 
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0.005 

0.007 

0.018 


0.075 

0.6 

0.002 

0.003 

0.005 

0.012 

0.018 
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Table 2. Stress intensity factors for internal edge cracks 
subjected to transient thermal stresses. (Icr(d) = 


z 

F 

IjgllggM 

R) = .0005 
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fb = .05 

kT(d) 
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0.117 


- 10 - 












Table 3. Stress Intensity factors for internal edge cracks 
subjected to transient thermal stresses, (kj(d)= 

■ ^ • F = 0-5- 9 .-T.-T 0 ) 


l 
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Fo=.0005 1 
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kj(d) 
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m 

m 


0.997 


M 

0.644 


mm 

0.364 

0.485 


mm 

0.525 


0.064 

0,146 

0.214 
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0.013 

m 
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Table 4. Stress intensity factors for internal edge cracks 
subjected to transient thermal stresses. (k,(d)= 

_ k(d) l-v a _ n T r-:_ Dt u_._ . . _t- t \ 

^ W » Qc-Tco-Tq) 


l 

h 

iSBsyi 
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0 

• 

II 
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Table 5. 


Stress intensity factors for internal edge cracks 
subjected to transient thermal stresses. CkjCd)“ 

• I" “•5- ^ «"T--To) 


F 

BSBiiiUB 


o’ 

II 

• 

o 

o 

Fo=.005 


kyCd) 

n<fCdT 

kyCd) 


WBM 

0.833 

WM 

0.277 

0.1 


0.701 


0,247 

0.2 

■IH 

0.589 
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0.501 
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Fig 2 Transient temperature distribution 0/9oo in a hollow 
cylinder due to a sudden temperature change on the 
inner radius. a/b=0,3, h=b-a» Fo=Dt /b^,0/0<» = 
(T(r,t)-To)/CT«-To). 
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(r-a)/h 


Fig. 3 Transient thermal stresses in a hollow cylinder which 
has been suddenly cooled by a temperature T« on its 
inner radius. a/b=0.3, h=b-a, Fo=Dt /b^, a* = 
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(r-a)/h 


4 Transient temperature distribution 9/0co in a hollow 
cylinder due to a sudden temperature change on the 
inner radius, a/b=0,5, h=b-a, Fo=Dt /b2,0/0co = 
(T(r,t)-To)/(T„-To). 
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(r-a)/h 

Fig. 6 Transient temperature distribution 0/0« in a hollow 
cylinder due to a sudden temperature change on the 
inner radius. a/b=0.7, h=b-a, Fo=Dt /b^.e/e^ = 
(T(r,t)-To)/(T„-To). 
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Transient thermal stresses in a hollow cylinder which 
has been suddenly cooled by a temperature T« on its 
inner radius, a/b = 0,7, h=b-a, Fo=Dt /b^, a* = 
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(r~a)/h 


Fig. 8 Transient tei^Jerature distribution e/e,, in a hollow 
cylinder due to a sudden temperature change on the 
inner radius. a/b=0.9, h=b-a, pQ=Dt /b^,0/0o9 = 
(T(r,t)-To)/(T«-To). 
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Transient thermal stresses in a hollow cylinder which 
has been suddenly cooled by a temperature T<» on its 
inner radius, a/b = 0.9, h=b-a, Fo=Dt /b^, o* = 
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